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ABSTRACT 

A method  to  manually  calculate  the  local  wind  current  in  deep  water 
is  developed.  The  method  is  based  on  time-dependent  Ekman  dynamics.  The 
method  is  tuned  and  tested  against  a 21/2  month  long  current  and  wind 


record  in  the  North  Atlantic. 
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INTRODUCTION 


The  local  wind  current  is  that  current  in  the  ocean  which  is  generated 
solely  by  the  action  of  the  local  wind  field  acting  on  the  ocean  surface. 

The  ability  to  calculate  this  local  wind  current  by  a manual  method  is  an 
important  aspect  in  SAR  Planning.  The  development  of  a manual  method 
presents  a two  fold  problem:  (1)  the  selection  of  an  appropriate  model 
to  determine  the  local  wind  current  and  (2)  the  adaption  of  the  model  to 
a manual  method. 

THE  MODEL 

Most  of  the  models  to  determine  the  local  wind  current  are  extensions  of 
Ekman's  (1905)  work.  Ekman's  (1905)  model  for  the  determination  of  the 
local  wind  current  for  an  unbounded  infinitely  deep  ocean  is  based  on  the 
following  assumptions: 

a.  The  fluid  is  hydrostatic 

b.  The  fluid  is  homogeneous 

c.  The  eddy  viscosity  is  constant 

d.  The  lateral  stresses  are  neglected 

e.  The  non-linear  interaction  terms  are  neglected 

With  these  assumptions  the  equations  of  motion  in  complex  notation  can  be 
written: 

where 

w = u+lv,  horizontal  components  of  velocity 
f * Coriolis  parameter 
V « vertical  eddy  viscosity  coefficient 
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This  equation  states  that  the  acceleration  of  the  fluid  is  balanced  by  the 
Coriolis  force  and  the  frictional  force  due  co  vertical  shearing  stresses. 

With  the  boundary  conditions  that  there  is  no  flow  at  the  bottom  of 
the  ocean,  that  the  stress  at  the  fluid's  surface  matches  the  wind  stress 
and  that  the  velocity  at  t^O  equals  0,  equation  Cl)  can  be  solved  by  the 
method  of  Laplace  transforms.  This  solution  as  derived  bv  JelesnlansVl 
(1970)  is  then: 

c 

where 

p^=  density  of  the  water 

H = depth  of  the  ocean 

F = Fx  + iFy,  complex  wind  stress 

Jelesnlanski's  (1970)  solution  was  selected  because  it  is  an  exact  solution 
to  the  equations  of  motion  under  the  stated  assumptions. 

To  Implement  equation  (2)  it  is  necessary  to  choose  a relationship 
between  the  wind  stress  and  the  wind  velocity.  For  ease  of  computation  it 
was  assumed  that  the  wind  stress  was  linearly  proportional  to  the  wind 
velocity,  i.e. 


F(T)  - PaCijWd) 

where 

Pg  ■ density  of  air 

Cq  ~ wind  drag  coefficient 

W ■ Wj5  + IWy,  complex  wind  velocity 
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Substituting  equation  (3)  into  equation  (2)^  a relation  between  the  local 
wind  current  and  the  local  wind  velocity  is  obtained: 

o 

Equation  (4)  essentially  states  that  the  present  wind  current  is  determined 
by  a weighted  average  of  the  previous  wind  history. 

ADAPTATION  OF  THE  MODEL  TO  A MANUAL  METHOD 
To  determine  the  local  wind  current, it  is  necessary  to  compute  the 
integral  in  equation  (4) . To  facilitate  this  computation,  assume  that  the 
wind  record  is  composed  of  a series  of  step  functions,  that  is,  the  wind  is 
assumed  constant  in  speed  and  direction  over  an  interval  2At  centered  at 
some  time  Tj . Equation  (4)  can  then  be  written: 


where 

N = t , number  of  time  intervals  in  the  wind  record. 

2A£ 

Now  write  the  complex  local  wind  current  as  a sum  of  a series  of  the  form: 


Substitution  of  equation  (6)  into  equation  (5)  and  elimination  of  the 
summation  signs  yields: 
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(7) 


Now  orient  the  axes  so  that  for  the  time  period  centered  at  Tj  the 
wind  blows  in  the  x direction.  Also  write  the  x and  y components  of  the 
complex  local  wind  current  contribution  for  the  time  period  center  at  Tj  as: 


Note  that  these  components  are  relative  to  the  direction  of  the  wind  over 
each  particular  time  interval.  Substitution  of  equations  (8)  and  (9)  into 
equation  (7)  yields  the  following  expressions  for  Cj and  Cj : 


(10) 


(11) 
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The  derivations  for  equations  (12)  and  (13)  are  given  in  Appendix  I. 

The  magnitude  of  the  local  wind  current  contribution  for  each  interval 


9 


Mfe) 


where 


If  aj  is  the  direction  from  which  the  wind  blew  over  the  interval  centered 


at  T j , then  the  direction  of  the  local  wind  current  contribution  for  that 
period  is: 
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(16) 


= oc^HS'cr-TA^-'^ 

Thus  from  a knowledge  of  C j (1)  and  Cj  (2)  the  contributions  to  the  local 
wind  current  can  be  determined. 

The  local  wind  current  can  be  determined  by  summing  up  the  contributions 
from  each  Interval: 

The  speed  and  Erection  of  the  local  wind  current  is  then: 


TESTING  AND  TUNING  OF  THE  METHOD 

Except  for  C^,  the  wind  drag  coefficient,  and  v,  the  vertical  eddy 
viscosity*  the  parameters  necessary  for  the  determinations  of  Kj  and  ijij  are 
accurately  known.  The  drag  coefficient  and  the  eddy  viscosity  can  thus  be 
used  to  tune  the  model. 

The  model  was  tuned  against  a time  series  of  wind  and  current  measure- 
ments made  with  current  meters  and  wind  recorders  by  the  Woods  Hole  Oceano- 
graphic Institution  during  June,  July  and  August  of  1970  in  the  Northwestern 
Atlantic  at  39°  07.6'N  and  70°  02.3"W  (see  Pollard  and  Tarbell,  1975).  The 
current  measurement  "ere  made  at  a depth  of  12  meters  in  water  of  depth  2682 
at  a site  designated  339D.  The  length  of  the  records  are  50  days,  20  hours, 
30  minutes. 
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The  wind  and  current  measurement  '-ere  averaged  over  6 hour  time 


periods  as  this  was  assumed  to  be  the  most  realistic  time  interval  for  which 
wind  measurements  could  be  obtained  in  a practical  situation.  Using  a time 
record  of  length  48  hours,  t = 48  hours  in  equations  (12)  and  (13),  Kj  and 
4ij  were  calculated  for  various  values  of  Cjj  and  v.  A time  record  of  48 
hours  makes  N = S.  The  various  Kjs  and  4ij's  were  then  used  in  equations  (17) 
through  (18)  to  reproduce  the  current  record.  The  distance  between  the 
calculated  drift  distance  and  the  observed  drift  distance  for  each  6 hour 
time  interval  was  then  used  as  a measure  of  the  accuracy  of  the  calculated 
velocities.  This  distance  can  be  calculated  from  the  equation: 

(21, 

where  U and  <{i  are  the  calculated  values  of  the  current  speed  and  direction 
and  Uq  and  <(1^  are  the  observed  values  of  the  current  speed  and  direction. 

A mlnmum  average  value  for  d was  obtained  for  Cp  = 1.9  x 10“^  cm/sec  and 
V = 50  cm2/sec.  These  values  are  realistic  values  for  the  drag  coefficient 
and  eddy  viscosity.  The  average  drift  error  for  these  values  was  1.5  nm. 

FORMUALTION  OF  THE  METHOD 

The  values  of  Kj  and  tj  car.  now  be  calculated  as  a function  of  latitude. 
These  calculated  values  are  shown  in  Table  1.  To  calculate  the  local  wind 
current  at  a specific  latitude  one  now  simply  has  to  multiply  the  local 
wind  speed  at  that  latitude  by  the  appropriate  Kj's  for  each  period  and  add 
the  appropriate  <)>j's  to  the  local  wind  direction  for  each  period.  The 
resultant  contributions  for  each  period  are  then  added  vectorially  to  obtain 
the  calculated  local  wind  current.  An  example  of  how  this  method  can  be 
used  is  given  in  Appendix  II. 
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APPENDIX  I 


Derivation  of  explicit  expressions  for  and 

The  expression  for  Cj from  equation  (10)  in  the  text  is: 


1-6 


The  expression  for  CJ 


can  be  obtained  by  taking  the  imaginary  part  of 


Equations  (1-6)  and  (1-7)  are  the  same  as  equations  (12)  and  (13)  in  the 


APPENDIX  II 


The  following  problem  Is  Intended  to  illustrate  that  various  computa- 
tional procedures  described  in  the  text. 

Given:  Present  time  is  2300Z  on  the  15th.  Information  is  received  that 
a distress  occurred  at  16008  on  the  15th.  Position  of  the  distress  is 
44®  -15'N,  58®-25'W.  The  search  is  expected  to  commence  at  0800Z  on  the 
16th  and  to  be  completed  by  16008  on  the  same  day.  The  following  wind 
information  is  obtained: 


Day 

jt 

Time 

Wind  Direction 

Wind  Speed 

16 

18008 

(“T) 

(kts) 

16 

12008 

340 

15 

16 

06008 

330 

20 

16 

00008 

320 

35 

16 

00008 

320 

30 

15 

18008 

260 

35 

15 

12008 

240 

30 

15 

06008 

230 

25 

15 

00008 

230 

20 

14 

18008 

230 

20 

14 

12008 

230 

20 

14 

06008 

220 

15 

14 

00008 

220 

15 
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In  order  to  calculate  the  local  wind  current  which  is  valid  for  the 
time  period  from  12008  to  18008  on  the  15th,  the  time  history  extending 
from  00008  on  the  14th  to  18008  on  the  15th  is  used.  The  appropriate  winds 
are  recorded  on  a wind  current  computation  sheet,  a typical  example  of  which 
is  shown  in  Figure  1.  The  column  labeled  45°N  latitude  is  then  located  in 
Table  1.  The  figures  from  this  column  are  recorded  on  the  wind  current 
computation  sheet  as  shown  in  Figure  1.  The  necessary  additions  and 
multiplications  are  then  made  to  determine  the  wind  current  contributions 
from  each  time  period.  For  example  the  wind  current  direction  appropriate 
for  period  1 is  obtained  by  the  following  addition: 

2e0°+22V  = 481“  = 121°T 

The  wind  current  speed  appropriate  for  period  1 is  obtained  by  the  following 
multiplication: 

35x0.023  = 0.80  knots 

These  values  are  recorded  on  the  wind  current  compulation  sheet.  For 
period  2 the  wind  current  direction  is  obtained  from  the  addition: 

240°+007°  = 247“! 

and  the  wind  current  speed  from  the  multiplication: 

30x0.010  = 0.30  knots 

Theve  values  are  recorded  on  wind  current  computation  sheet.  This  pro- 
cedure is  repeated  for  periods  3 through  8: 

Period  3 WC  direction  = 230“+136“  = 366“  = 006“T 

WC  speed  = 20x0,007  = 0.14  knots 

Period  4 WC  direction  = 230“+264“  = 494“  = 134“T 

WC  speed  =■  20x0.006  = 0.12  knots 

Period  5 WC  direction  = 230“+031“  = 261“T 

WC  speed  » 20x0.005  “ 0.10  knots 
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Period  6 WC  direction  - 230‘’+150°  = 389®  = 029“T 

WC  speed  = 20x0.005  = 0.10  knots 

Period  7 WC  direction  = 220°+286°  = 506°  = 146°T 

WC  speed  = 15x0.004  = 0.06  knots 

Period  8 WC  direction  = 220°+053°  = 273°T 

WC  speed  = 15x0.004  = 0.06  knots 

These  values  are  recorded  on  the  wind  current  computations  sheet.  There  is 
now  a wind  current  direction  and  speed  for  all  eight  periods.  These  values 
must  be  added  vectorlally  in  order  to  obtain  the  wind  current  direction  and 
speed  which  is  valid  during  period  1.  It  should  be  noted  that  the  wind 
current  direction  and  speed  for  each  period  is  not  the  wind  current  which 
existed  during  that  period,  but  is  only  the  contribution  from  that  time 
period  to  the  wind  current  which  is  valid  during  period  1. 

There  are  two  possible  methods  to  add  the  contributions.  The  first 
method  is  to  add  the  vectors  head  to  tail  on  a maneuvering  board.  The 
second  method,  which  is  convenient  if  a calculator  is  available,  is  to 
separate  each  contribution  into  its  east-west  and  north-south  components, 
add  these  components  and  then  reassemble  the  components  into  a vector.  In 
both  cases  the  resultant  vector  is  the  local  wind  current  valid  for  period 
1.  Both  methods  yield  a resultant  vector  of  136°  at  0.63  knots. 

In  order  to  calculate  the  local  wind  current  which  is  valid  for 
the  time  period  from  1800Z  on  the  15th  to  OOOOZ  on  the  16th,  the  time 
history  extending  from  06008  on  the  14th  to  00008  on  the  16th  is  used. 
Employing  procedures  the  same  as  those  described  above  a local  wind  current 
of  203°  at  0.67  knots  is  calculated. 
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similar  procedures  can  be  used  to  calculate  the  local  wind  currents 
for  the  time  periods  OOOOS  to  0600S  on  the  16th,  06008  to  12008  on  the 


16th  and  12008  to  18008  on  the  16th.  The  results  of  the  calculations  for 
all  the  desired  time  periods  are: 


Day 

Time  Period 

Current  Direction 

Current  Speed 

(“T) 

(kts) 

16 

12008-18008 

228 

0.27 

16 

06008-12008 

204 

0.36 

16 

00008-06008 

185 

0.49 

15 

18008-00008 

203 

0.67 

15 

12008-18008 

136 

0.63 

This  Is 

all  the  data  necessary 

to  calculate  the  drift  due 

to  the  local  wind 

current 

during  the  period  of  Interest. 
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Flaure  1 Wind  Current  Computation  Sheet 
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Table  1.  Coefficients  relating  drift  current  speed 
and  direction  to  wind  speed  and  direction. 
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TABLE  1 (cont) 


